CLASSICAL INTEGRABLE FIELD THEORIES IN DISCRETE 2 + 1 
DIMENSIONAL SPACE-TIME. 



SERGEY M. SERGEEV 

Abstract. We study "circular net" (discrete orthogonal net) equations for angular data 
generalized by external spectral parameters. A criterion defining physical regimes of these 
Hamiltonian equations is the reality of Lagrangian density. There are four distinct regimes 
for fields and spectral parameters classified by four types of spherical or hyperbolic triangles. 
Non-zero external spectral parameters provide the existence of field-theoretical ground states 
and soliton excitations. Spectral parameters of a spherical triangle correspond to a statistical 
mechanics; spectral parameters of hyperbolic triangles correspond to three different field 
theories with massless anisotropic dispersion relations. 



Introduction 

The "circular-" or "conic net" (or discrete orthogonal net) equations for angular data [4,5,8] 
take a selected place among all the classical integrable systems [3] on cubic lattice with 
AKP-type hierarchy. Algebraically, these equations arise as a Hamiltonian form of discrete 
three- wave system [6,20]. The "conic net" equations are classical g — > 1 limit of quantum 
"g-oscillator" model [2] - the top of a pyramid of three-dimensional quantum models - what 
guarantees in classics the existence of Lagrangian density, energy /action and variational prin- 
ciple [1]. The existence of quantum counterpart is an evident advantage of discrete space-time 
models with respect to their continuous space-time predecessors [22,23]. A straightforward 
geometrical condition for conic net equations is the reality of angular dynamical variables [1] 
of a circular net in Euclidean target space or of an ortho-chronous hyperbolic net in Minkowski 
onj^. However, this discrete differential geometry conditions can be essentially extended by a 
"physical" condition of reality of action near equilibrium point. 

A general complex solution of any AKP-type system in finite volume or a solution of Cauchy 
problem with generic initial data involves the algebraic geometry [15]. Such general solution 
of discrete "generalized conic net" equations in finite volume is known for a long time [9-14]. 
It involves a flat algebraic spectral curve Tg of genus g < {N — 1)^ for a size N^^ cubic lattice 
(three-periodic boundary conditions), ©-functions on Jacobian of Tg, and spectral parameters 
- three meromorphic functions on F^. A reduction of Tg to a sphere gives a g-soliton solution 
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(plane wave solitons) [17,18]. The soliton regime is the field-theoretical one since solitons are 
continuous excitations over a ground state - zero soliton homogeneous solution of equations of 
motion. Spectral parameters in the soliton regime are a triple of complex numbers, they enter 
directly into equations of motion as extra parameters providing the existence of homogeneous 
solution. The spectral parameters break the straightforward discrete-geometric interpretation 
of equations of motion. 

These two principles - the reality of energy/action and existence of homogeneous solution 
for "non-zero" spectral parameters - are starting points for classification of physical field the- 
ories and statistical mechanics for "generalized conic net" equations. There are four distinct 
regimes of spectral parameters and corresponding regimes of dynamic fields providing the 
reality of action. Parameterizations of ground states have a structure of cosine theorems for 
spherical or various hyperbolic triangles. In the case of spherical triangle the ground state is 
the absolute minimum of energy functional and thus it corresponds to statistical mechanics. 
In three field-theoretical cases of hyperbolic triangles a solution of equations of motion pro- 
vides an extremum of action; a value of whole action on soliton solution does not depend on 
amplitudes of solitons. Expressions for soliton waves involve projective coordinates of hyper- 
bolic triangles, a general field-theoretical solution of equations of motion is a set of soliton-anti 
soliton pairs analogues to elementary stationary waves. In low energy-momentum limit the 
soliton plane waves have a cone-type (anisotropic massless) dispersion relation. 

All the regimes of spectral parameters and dynamical fields have manifest counterparts in 
quantum case. Regimes of fields correspond to classical limits of different representations of 
g-oscillator algebra; related regimes of spectral parameters correspond to either real or unitary 
quantum i?-matrices. However, relations between quantum theories and classical theories are 
not straightforward. 

This paper is organized as follows. In Section 1 we fix notations for the "conic net" equations 
generalized by spectral parameters. Following [1], Lagrangian density is defined in Section 
2. Next, in Section 3 we classify ground states. Soliton solutions of equations of motion and 
dispersion relations are defined in Section 4. In Section 5 we describe quantum counterparts 
of our four regimes. Finally, in concluding Section 6 we discuss roughly a place of finite gap 
solutions. 



1. Generalized conic net equations. 



Let n be a node of a big size simple cubic lattice 



(1) 



n = (ni, 7x2,7x3) , riieZ . 



Let ei, 62, 63 be the unit vectors for the lattice 



(2) 



ei = (1,0,0), 62 = (0,1,0), 63 = (0,0,1), 
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SO that n = niCi + n2e2 + ^.363. With each (n — ei,n) edge of the cubic lattice we associate 
a doublet of dynamical variables 

(3) Ai,n = ik^n^ af^J , such that fc^^ =^ 1 - a+^ar^ . 

Last relation here is in fact the "conic net" condition. The local equations of motion relate 
the neighbors of every node n, 

(4) Ai = Ai^n on (n-ei,n) and A'i = Ai^n+Si on (n,n + ej), 
as follows: 

(5) {afy = afaf - u^^k^k^af , 

ik2afy = uf^ {k^af + u'^^k^a'faf) , 

and 

(6) ^1^2 — ^1^2 1 ^2^3 — ^2^3 ■ 

C-valued parameters Uj - exponents of spectral parameters - are the same for all nodes n. 

This classical system can be viewed as an extension of discrete three-wave equations and 
conic nets since the last ones correspond to trivial spectral parameters. The circular net 
in Euclidean geometry is described by regime k"^ > 0, real and ui = —U2 = = 1. 
The Euclidean "circular net" point has a smooth continuous limit ~ the classical three-wave 
resonant equations [23]. The hyperbolic net in Minkowski geometry is described by regime 
/c^ < 0, real and ui = U2 = = 1- For non-trivial spectral parameters or for complex 
a geometrical interpretation of equations ^ in unclear. 

However, relaxing the geometric condition, equations ([5]) define an evolution in 2 + 1 di- 
mensional space-time. Discrete time is t = ni + 71-2 + n3, so that ([5]) literally gives the map 
from time t to time t + 1. A straightforward way to introduce space-like coordinates is to take 
rii = X and = y so that 

(7) n = ni(ei - 62) + 723(63 - 62) + (ni -I- 71-2 + 713)62 . 

The Cauchy problem is well posed for finite size space-like surface, 

(8) x,yGZAr, iV>l, tGZ>o. 



The evolution corresponds to a relativistic field theory since the locality of evolution map pro- 
vides the relativistic casuality. Note, there is no way to introduce a usual local Hamiltonian 
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for classical discrete time evolution; for instance a principal object of corresponding quan- 
tum theories is a discrete time Heisenberg evolution operator. The framework of statistical 
mechanics implies 3D Dirichlet or 3-D periodical boundary conditions for system ([5|), 

(9) niGZ^v, i = 1,2,3, iV>l. 

The principal difference between a statistical mechanics and a field theory is that for given 
reality regime and 3D periodical boundary conditions a statistical mechanics has an unique 
ground state minimizing an energy while a field theory with a saddle-type action has a rich 
structure of stationary modes. 

Equations ([5]) is a canonical transformation preserving locally the g-oscillator symplectic 
form [2], 

^^Q^ daf Ada- ^ ^ (ia+ Ada-" 

i=i i=l ^'i 

and thus they have the discrete-type Hamiltonian structure [1]: 

(11) ^oz\K\ = 2v[^G{v]v') , \og\ki\ = --Vi—G{v]v') , 

where 

def af 



(12) 



is a useful canonical partner to ki, the variables Vj,v'j, kj, k'j are related by ([5|), and G{v;v') 
is a generating function of the map ^ . 

It is more convenient to treat fe, k' related by ([6]) as independent variables and define G{k;k') 

by 

3 

2 



1 ^ 

(13) dG{k-k') = -Y,{^og[v[]d\ogk[-\og[v,]d\ogh) . 



i=l 

Here 



(14) [v] = V if t> is positive or unitary and [v] = —v if v is real negative. 

Negative v corresponds to regime k'^ > 1. Functions G{k;k') and G{v;v') are related by the 
Legendre transform, 

(15) G{v-v') = \ {^\og\v'^,\og%\-\og\v,\\og\k^}\^ - G{k;k') . 

Function G{k; k') is preferable since ki are the field coordinate-type variables with fixed reality 
regime A;^ G M. 
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The sum of local generating functions over all nodes, 

(16) si ~ G{Vi^n\Vi^n+e.^ , 

gives an action/energy for whole lattice. Equations of motion are the extremum conditions 
for (dH) [1]. 

Considering the equations of motion as the quasiclassical limit of quantum models, we 
expect two distinct regimes: A:'^ < for modular representation of g-oscillators and /c^ > for 
Fock space representations (regime of unitary k for cyclic representations we do not consider 
here). Also, it is evident from (jl3p that si has two regimes of reality: either the regime of 
real v when G is manifestly real or the regime of unitary v when iG is real. Thus we expect 
a "priori the existence of four distinct regimes: 

Regime 1: /cf < , \vi\ = 1 , Regime 2: k^KO, G R , 
Regime 3: k^ > , \vi\ = 1 , Regime 4: k^ > , f , G M . 

The presence of external parameters Ui in equations of motion is of great importance for 
classification scheme because of the regime f j E M corresponds to G M and the regime 
\vi\ = 1 corresponds to \ui\ = 1. 

Note also a signs symmetry of equations ([5]): a change of sign of any two parameters of ut 
followed by a change of signs of corresponding af^ on e.g. even edges (n^ in n = (ni,n2,n3) 
- even) does not change the whole set of equations of motion. This transformation preserves 
the variables Vi^n- Contrary to this, a change of a sign of a single Ui changes a structure of 
solution of equations of motion substantially. 



2. Generating function 
The {v, A;)-form of local equations of motion ^ is the following [1]. Let 



^0 — / / ^1^3 5 ^1 — / 2 ' 2 — ; ^3 — / 2 ' 

, _ V2 I , _ v'^ 2 , _ ujul , _ v'2 „ 

''"0 — 2 ' '^i ~ T % ' '^2 — ~r~i 2~ ' ''"3 ~ ~ ""i ' 



;i8) 



^1^3 "2 ''l''3 ''1''3 "2 ""l^S 

so that TqTq = t^t[ and 

(19) TjT- = t'j.t'i , i, j, fc, I = any permutation of (0, 1, 2, 3) 
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Equations ([5]) provide 



(20) 



^3 



(1 


-T2)(l 


-^3) 


fcf = 


(1 




-^3) 


(1 






(1 


-^o)(l 


-r^) 


(1 






^2^ = 


(1 


-^o)(l 


-^2) 


(1 


-Ti)(l 


--3) 


(1 


-r[){l 


--^) 


(1 


-Ti)(l 






(1 






(1 






(1 


-^o)(l 





Note, expressions for kf^k'^ are invariant with respect to inversion r^- 



Thus, the differential (jl3p has the form 

3 

4 



1 ^ 

A;') = t ^ (log[r,]dlog(l - r,) - log[rj]dlog(l - rj)) 



(21) 



+ 2 (loguidlog/c'i - logn2dlogA;2 + log Ugdlog ^3) . 
Integrating this, one gets 

(22) G{k- k') = Go{k; k') + ^ (logn? log \k[\ - logu^ bg 1^1+ logu^ log , 
where 

Go{hk') = -Y.{J[r^)-J{r'^)) 



(23) 



and 



(24) 



i=o 



J(r) 



20 



log 2; (ilog(l — z) in Regimes 1,2 

■r 

logz (ilog(l + z) in Regimes 3,4. 



20 



Choice of zq common for all integrals is irrelevant. We choose zero value for a constant of 
integration in ()22p . The GQ{k]k')-iexva. in (j22p corresponds to Lagrangian density without 
fields [1] but it must be treated carefully due to inversion symmetry of Tj. 

Generation function (j22p becomes symmetric after an elementary gauge transformation, 

(25) Gsymik- k') = Go{k; k') + ^ (logu? log\k^k[\ - logul log {k^k'^l + loguj log {k^k'^^l) . 
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3. Ground states 
By "ground state" we understand the homogeneous solutions of equations ([5]): 

(26) = v] , = k'j , 
what corresponds to 

(27) J^ = J^^r; = rrK 

In the vicinity of this point the fields are parameterized by 

^2 ^ (nin3 - U2){uiU2 - U3) ^2p, jj2 ^ (-^1^3 - U2){uiU2 - U3) ^^p', 
^ -■ -• ^ '1 {1 - UIU2U3){U2U3 - Ui) 

-,2 _ \^ ~ "-r"-2"-3n"-r"'3 ~ "-2; 2p2 i,'2 — (1 ~ '"1^2^3)(^1^3 " ^2) 2p^ 



(28) ki = 7 -i-z-.M-i-o ^z, 2p, ^ 



(nin3 - U2){uiU2 - 


- U3) 




(1 - liiU2U3)(w2li3 


- Ul) 


1 


(1 - UiU2U3){uiU3 


- U2] 




{U2U3 - Ui){uiU2 - 


- U3) 




{U2U3 - Ui){uiU3 - 


- U2) 




(1 - UiU2U3){uiU2 


- U3) 




Pi 


-A 


= (^2 



^2 ^ {U2U3 - Ui){uiU3 - U2) 2p, jj2 ^ (^2^3 - Ui){uiU3 - U2) Op', 



(U2n3 - Ui){uiU2 - U3) 

{U2U3 - Ui){uiU3 - U2) 
(1 - UiU2U3){uiU2 - U3) 



where 
(29) 

The generating function (|25p near the ground state /? = is 



Gsym — Gsy'm\p=p'=0 ~l~ ^ ^ 



„2 



(30) 



lil + til ^ Uo + Uo^ Un+ tin ^ 

_i "''2-'^3 ~r -1 "'^1'''3 ~ -1 •^1-^2 

lii — til Ur, — Un tiq — tiq 



(l(lU2ti3 - l)(ttlti2 - n3)('Wl^i3 - ^i2)(^i2^t3 - Itl) .2 , ^/^^3^ 



where 
(31) 



ufululiU;^ - «i ^)(n2 - U2 ^){U3 - U3 ^) 



Pi +Pi r2 - Pi 



5' + 0{p' 



2 



and Gsym\p=p'=o = G{k; k) is calculated at the point = p'^ = {k^ = k[). 

Now we are ready to classify the generating functions for all four regimes. Depending on 
regime, expressions ()28p are equivalent to cosine theorems for spherical or hyperbolic triangles; 
thus the classification scheme is based on the spherical and hyperbolic geometry. 
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3.1. Regime 1. Unitary spectral parameters are given by 

(32) ui = e'^K^i , U2 = e''2</'2 , = e'^3</'3 ^ 
where (^j > and are signs. In what follows we use short notations 

(33) /34 = ^ , /3j = /34 - . 

In spherical geometry the angles are sides of a spherical triangle. It is convenient to define 
excess /5o instead of half-perimeter /34 by 

(34) /Jo = vr - /?4 . 
Then for ei = €2 = €3 = 1 one has 

(35) ^1 = - tan^ y , fc^ = - cot^ y , = - tan^ y , 

where di are dihedral angles of the spherical triangle with sides : 

cos 4>i — cos 4>j cos </)fc cos di + cos 0,- cos 

cos = . , . , , cos 4>i = ■ Q ■ a ' 

sm (pj sm (pk sm Oj sm 

(36) 

2 sin/3, sin/3fc 

tan — = — . 

2 sm po sm pj 

Spherical geometry implies positive ^-excess 

(37) ei + 92 + 9z>TT => 0<Pj<TT, j = 0,1,2,3. 

Other choices of signs such that €162^3 = 1 are equivalent to crossing transformations of 
spherical triangle (in fact, crossing transformations involve (/> — > vr — the sign symmetry of 
Ui we discussed above). 
For arbitrary signs Ci let 



(38) F{k;k') = \eie2e3Gik;k') , JT = ^ ^(A;^,^; 



Homogeneous solution provides the absolute minimum of functional J^; on this ground state 
the free energy density F{k; k) = Fq is given by 

3 3 

(39) Po = Y. -"(^i) = Yl - -^(^4) > , 

i=o i=l 

where the Milnor's Lobachevski function is 

(40) Ji(p) = - log(2sinx)dx . 

Jo 

The statement about absolute minimum can be verified instantly in the free field approxima- 
tion (l30l) where 



(41) = N^Fq + positively defined quadratic form of pi 
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Here A^^ is a volume of the lattice. We will discuss this statement beyond the free field 
approximation in the next section. 

3.2. Regime 2. This is the case of spectral parameters 

(42) ui = e^i'^i , U2 = e^^"^^ , U3 = e^^'^^ , 

where (pi > and are signs. Values of ki of homogeneous solution for ei = £2 = £3 = 1 
are given by (j35p where 9i are dihedral angles of a triangle on upper sheet of two-sheets 
hyperboloid. Parameters (pi are hyperbolic sides of this triangle. Cosine theorems read 

cosh (pj cosh — cosh (pi cos 6i + cos 9j cos 9k 

= — ^"n^ ' = — ^ ' 

smh (pj smh (p^ sm Oj sm t^fe 

(43) 

n 0j sinh /?,■ sinh (3k 

tan — = —r — — - , 

2 sinh p4 sinh pj 

where excesses are given by (I33p . Hyperbolic geometry implies negative 0-excess 

(44) 9i + d2 + es<TT 0<ft,z = l,2,3. 

Other choices of signs are analogues of crossing transformation of the hyperbolic triangle. 
For arbitrary signs define the Lagrangian density and the action by 

(45) L{k;k') = -eie2e3G{k;k') , = L{ki^n; h^n+e,) ■ 

The criterion for a correct sign of Lagrangian density is the positive sign near d'^ in the 
free-field approximation ([30]) . 

(46) L{k- k') = (positive coeff.) x 5^ - V{x) - Vq , 

so that 5^ stands for a square of velocity and V{x) + Vq stands for a potential. In this regime 
the quadratic form V{x) is positively defined. On homogeneous solution (ground state) 

(47) L{k;k) 



"^0= / log(2 sinh x) dx —y / log(2 sinh x) > 0. 
^0 Jo 

The global quadratic form has saddle structure, general solution of linearized equations of 
motion are plane waves with a certain dispersion relation. As it is clear for free theory, the 
value of the whole action on any plane wave solution of equations of motion in finite volume 
coincides with its value on vacuum solution, 

(48) A = -N^ Vq . 

In the next sections we give a general solitonic solution of field-theoretical equations of 
motion which can be regarded as excitations over the ground state. Dispersion relation for 
solitons is the same as the dispersion relation for linearized theory. Since for 3D periodical 
boundary conditions the value of whole action at equilibrium point is an universal invariant. 
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it does not depend on amplitudes of solitons and therefore it equals to value of the whole 
action for ground state. 

3.3. Regime 3. Unitary parameters Ui are given by 

(49) Ml = e'^i'^i , U2 = -e'^2<^2 ^ ^ gi.3</'3 

where as usual > 0, ei are the signs and the excesses are defined by (p3|) . This gives for 
ei = £2 = £3 = 1 

(50) kl = coth^ ^ , kj = tanh=^ y , kj = coth^ y , 
with cosine theorems 

cos + cos cos (pk cosh cosh 6*^ — cosh 6i 

cosh^i = : ; r^^- , cos (Pi = - 



(51) 



sinh 9j sinh 9^ 



coth^ - 



2 COS Pj COS 



2 cos Pi cos /3j 

This is a hyperbolic triangle formed by an intersection of three planes with time-like normals 
(and hyperbolic angles 6i between them) and one-sheet hyperboloid. Trigonometric sides (pi 
are defined in motionless frame of reference for each plane. The time-like normals form a dual 
triangle on two-sheets hyperboloid of Regime 2. The geometry provides the constraint for 

(52) 0</3i<7r/2, i = 1,2,3,4, 

otherwise it would be Regime 1. Other choices of signs ej are analogues of crossing transfor- 
mation of the hyperbolic triangle. 

We define the Lagrangian density for arbitrary signs by 

(53) Lik;k') = kie2€3Gik;k') , 

where the sign criterion is the same as for Regime 2. However, the quadratic potential here 
is not sign defined. On homogeneous solution (ground state) 

(54) L{k;k) = -Vo= / log(2 cos x) - V] / log(2 cos x) < 0. 

Jo Jo 

3.4. Regime 4. Real spectral parameters are defined by 

(55) m = , U2 = ^ ^ ^-es^s ^ 

where (pi are positive, are again signs. Negative sign of one of Ui makes the difference with 
Regime 2. This parameterization gives for ei = e2 = £3 = 1 

(56) A:2=tanh2^, A;| = coth^ ^ , fef = tanh^ ^ 
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where cosine theorems are 

cosh ipi + cosh (j)j cosh (j)]. cosh Oi + cosh 9j cosh 9^ 

cosh Oi = , cosh (t)i = . I , 

smh (pj smh ip^ smh 9j smh t^fc 

(57) 

^ ^2 cosh /?j cosh 

2 cosh cosh /3j 

The hyperbohc excesses are defined by (|33|) . This corresponds to a triangle on one-sheet 
hyperboloid with hyperbohc sides 0j and hyperbohc dihedral angles Such triangle is 
the section of one-sheet hyperboloid by planes with space-like normals. Note, two planes 
with space-like normals and hyperbolic angles between them do not intersect on two-sheet 
hyperboloid. Contrary, two planes with space-like normals and trigonometric angle between 
them intersect on two-sheets hyperboloid, this corresponds to Regime 2. 

Other choices of signs are analogues of crossing transformation of the hyperbolic triangle. 

Lagrangian density and the action are then 

(58) L(/c;A;') = eie2e3G(/c;A;') , .s^= ^ L(/ei,^; A;i,„+eJ . 

In this regime the quadratic potential is not sign defined either. Ground state gives 

./34 3 

(59) L{k\k) = -Vq = / log(2coshj;) dx - > / log(2 coshx) >0. 

h i=i-^° 

4. SOLITONS 

4.1. General soliton solution of ([5]). A general (complex) soliton solution of equations (l5|) 
is given by the reduction of general algebraic geometry solution corresponding to the reduction 
of genus g curve to a sphere with punches. Resulting expressions are the following [17]. For 
the number of solitons 5 > let 

(60) {^i,^.}i=l,...,9, / = {/.}i=1..5 

be a set of 8(7 complex values. For given / and {Xj^Yj) let 
Define next 

aei \j\ 

(62) G(/) 



W^>,{X^-X,) 

where in the numerator there is the determinant of g x g matrix with matrix indices j, k. In 
this expression fj is the amplitude of jth soliton, if one of fj = then (j62p simply gives g — 1 
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soliton expression. For instance, 

5 = ^ e = i, 



(63) 5 = 1 => e(/i) = i + /i, 

o ^ mff f\ -1 , f , f , (Xi - X2)(Yi - I2) 

5 = 2 ^ e(/l, /2) = 1 + /l + /2 + _ y^^^y^ _ , 

etc. In general, at the first order of / 

9 

(64) e(/) = 1 + ^ + higher terms, 

what corresponds to the free field (linear) approximation. Let further 

(65) = MXj,Y,) = I^.I p^^|(y?' I Q j . j = h-,9, k = l,2,3 



and 



(66) /(n) = {/,(n)},=i„„,,, /,.(n) = c^^W^^c^g"^ 



Let also for brevity 

(67) G„ = e(/(n)) . 

General soliton solution of ([5]) is then given by [13] 

-62+63 



(68) 



Tl— 62 +63 

2 _ E{Qi,Q3)E{Pi,P3) @r.~., 0TI+61-62+63 



2,Tl 



-E'(Ql5 -P3)-£'(-Pl) Qs) 0T1+61 -62 ©71-62+63 
2 E{Ql,Q2)E{Pl, P2) ©T10TI+61-62 



'>'3,Tl 



where 

(69) E{Q,P) 



E{Qi,P2)E{Pi,Q2) e^+eiO 

Tl— 62 



./dQdP 

is the prime form on a compact complex plane. Spectral parameters in this parameterization 
are given by 

„ „ „ _ E{Ql,Q2)E{Pl,Q3)E{P2,Ps) Ml _ i;(Qi.P2)g(Pi,P3)i?(Q2,Q3) 

"1"2"3- E{Pi,P2)E{Qi,P3)E{Q2,Q3) ' «2«3~ E{PuQ2)E{Qi,Q3)E{P2,P3) ' 

(70) 

M2 = E{Pi,Q2)E{Qi,P3)E{P2,Q3) U3 ^ E{Pi,P2)EiQi,Q3)E{Q2,P3) 

U1U2 E{QuP2)EiPi,Q3)E{Q2,P3) ' «i«2 E{Qi,Q2)E{Pi,P3)EiP2,Q3) ■ 
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4.2. Identification of parameterizations. Homogeneous solution of ([5]) corresponds to 
g = when all Qn = 1 • Expressions (j68l) are equivalent to parameterizations of homogeneous 



kj in terms of spherical and hyperbolic triangles. Let us demonstrate this statement in more 
details. 

In Regime 1 of Euclidean spherical trigonometry, consider thee planes defined by their unit 
normal vectors Hi in some auxiliary frame of reference, 

( 71 ) fii = (sin -i? j cos ifi , sin ■!? j sin ipi , cos i?i ) 

The dihedral angle between two planes equals to the angle between normals, 

(72) 6*3 = nin2 , 02 = n - nin^ , Oi = n2n^ , 

where 9i are inner dihedral angles of the spherical triangle. Cosine theorems give 

(73) cos 6*1 = (n25 "^s) = cos '!92 cos -i^s + sin'!92 sin??3 cos(v32 — V's) ^tc. 

For all points fii on the sphere define their stereographic projections to a complex plane: 

(74) Qi = tan ^e'^' , Pi = - cot ^e'^» . 
The cosine theorem can be then rewritten as 

_. 2O1 ^ COS^i-l ^ {Q2 - Q3){P2 - P3) 

^ ' 2 cosOi + l {Q2-P?.){P2-Q^) ' 

what makes exact correspondence between (p5|) and ([68|) for 0„ = 1. 

This can be done similarly for all other regimes. In Minkowski metric g = diag(l, —1, —1), 
a time-like unit vector is parameterized by 

fii = (cosh -di , sinh -di cos ipi , sinh t?j sin ipi) , so that 

(76) 

cosh 6*1 = (n2i ^^3) = cosh ■!?! cosh t?2 — sinh -di sinh t?2 cos (991 — v?2) , 
what gives for Regime 3 

(77) Qi = tanh y e'^^ , Pi = coth y e''^' . 

The proper parameterization of space-like unit vectors is the following: 
fii = (tan i^j, sec 'di cos v^j, sec f}i sin (pi) , so that 

(78) 

a .0 ^\ cos((/?2 - 953) - sini?i Sin'!?2 
costal or coshWi = -(n2,n3) = , 

cos tT2 COS 173 

what gives for Regimes 2,4 

(79) Qi = e'(^'+''») , Pi = -e'^^-^^) . 
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Thus, in terms of complex parameters Pi,Qi, Regimes are classified as follows: 

Regime 1: P^Q* = -1 , 

(80) Regime 3: P^Q* = 1 , 

Regimes 2,4: \Pi\ = \Qi\ = I , 

where * stands for complex conjugation. 

4.3. Plane waves and dispersion relation. Relation (j66p stands for a plane wave with 
exponential frequencies uJk{X,Y), k = 1,2,3. Parameterization (|65p can be viewed as a 
general solution of an algebraic equation relating three w^, = 1, 2, 3. This dispersion relation 
can be obtained by elimination of X,Y from ()65p . it has the form 

2 

(81) ^ Cijk{ui,U2,u:i)uj\uj^2^l = {) , 

i,j,k=0 

where Cij^k{uiiU2iU3) are simple but lengthy rational coefficients. Note that the free field 
approximation provides the same dispersion relation. 

Three-periodical boundary conditions in a rather big volume require unitary uj^. Parame- 
terization (j65p and definition of regimes ()80p provide immediately the following: 

• In Regime 1 it is impossibl^ to make all three u;k{X,Y) unitary. Thus, the homo- 
geneous solution is indeed the absolute minimum of the energy functional (j38p . For 
the open boundary conditions, the solitons of Regime 1 break the signature condition 

< 0. 

• In Regimes 2,4, when Pk,Qk are unitary, all uJk{X,Y) are unitary if X*Y = 1. 

• In Regime 3, when P^-Ql- = 1, all u]k{X^ Y) are unitary if |X| = |y | = 1. 

In all field-theoretical regimes the reality condition = 0„ is satisfied for soliton-antisoliton 
pairs with conjugated amplitudes. 

The dispersion relation for near unity, 

(82) oji = e'P' , u;2 = 6"'^^ , 003 = e'P' , 
where momenta pi are small, reads 

(83) PiP^iU'i — ^)^ ~^ PiP2P\{'^i ~ ■'^r^)('"2 ~ ""2 "'^)(% + % "'^) + cyclic permutations = 0. 

Due to the homogeneouty, this relation describes a cone-type surface in momentum space. In 
the symmetric cases ui = U2 = U3 or ui = —U2 = U3, what corresponds to = (1)2 = (ps = 
(p > 0, relation ([83]) becomes 

(84) pIpI + pIpI + pIpI + 2CpiP2P3{pi +P2+P3) = , 



^Unitarity condition in Regime 1 demands X*X = Y*Y 
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where 

Regime 1: C = cose/), , -1/2 < C < 1 ; 

Regime 2: C = cosh0 , 1 < C ; 

Regimes: C = -cos0, -1< C < -1/2 ; 

Regime 4: C = — coshc/) , C < — 1 . 

Our four regimes cover the real axis, C G M\{1, — 1, — 1/2}. Define then an energy E and 
space-hke momenta tTj (we do not care about scales of energy and momenta) by 

(86) Pi = E + TTi , VTl + 7r2 + VTs = . 

Let 

(87) vr = -(vTi +7r2 +7r3) , 7 = ^3 ■ 

If C < — 1/2 or 1 < C, equation ([84|) defines an anisotropic cone- type surface 

(88) E = a(7)7r 
where 0(7) is a real solution of 

(89) {C + ha^-Ca^ + iC-l)ja + l = 0. 
Anisotropy parameter 7 is bounded, 

'4 



(90) - 70 < 7 < 70 , 70 w 27 
Critical values 7 = ±70 correspond to three selected directions in the momentum space when 

(91) p^=p2 = or = p3 = or p2 = = . 

When C < —1/2 (Regimes 3,4), equation has one positive and one negative solution 
what gives a rather anisotropic "cone" with 




where a+ij) are taken positive (negative solutions for given 7 are 0^(7) = — a+(— 7)). When 
7 = ±7o, equation ([89l) has extra solutions a = it-^, these solutions are isolated and therefore 
do not belong to one-parameter family, they have no relation to and should be ignored. 

When 1 < C (Regime 2), equation ()89p has two positive and two negative solutions what 
gives two imbedded tangent anisotropic "cones". The "cones" are tangent along 7 = ±70, 
a(ib7o) = i"^) these points are not isolated. Existence of two "speeds of lights" is a surprise. 

Regime 2 involves the Lorentz group limit. If all in this Regime are small, cos(^i ~ 1, 
then the dispersion relation ()83p becomes 

(93) (PiP2^PiPs^P2P^y ^ 
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what is equivalent to isotropic light cone and gives the pure Minkowski metric in the momen- 
tum space. 

When —1/2 < C < 1 (Regime 1), equation ([89|) has no real solutions as expected. 

5. Quantum theories 

In this section we discuss the relation of classical regimes to quantum models. 
We commence with a short remainder of quantum i2-matrices. Let A be the enveloping of 
the g-oscillator algebra 

(94) ka^ = q^^a^k , a+a = 1 - q'^k"^ , a = I - qk^ . 
equipped by a pair of C-valued parameters A, ^, 

(95) ^= (l,fe,a±,A,/i) . 

The map -R123 [1, 2] of tensor cube Ai® A2® A^ defined by (confer with ([S])) 

-^123 ^2^t -^123 ~ ^1^^ (^3*^1^ + '"2'^^l"2 *^3 ) ' 

(96) -^123 ^2 -^123 ~ *^3^ ~ ^2'^^1^3^f ' 
-^123 ^2^3 -^123 ~ ^3^^ (^1^3^ + '"2 ^^3^? *^2 ) ' 

where 

(97) m = — , U2 = , ns = — , 

A2 A1//3 /i2 

satisfies the adjoint tetrahedron equation in A®^ and the quantum tetrahedron equation in 
proper Rep(^)®^ (spectral parameters Ai, ^ui, . . . Ae, /^e foi' the tetrahedron equation are free). 
Equations (j96p provide in addition 



(98) -Ri23^1^2-f^l23 ~ ^1^2 5 -^123^2^3-'^123 ~ ^2^3 ■ 

"Constant" matrix ri23 corresponds to ui = U2 = = 1. In modular representation 

(99) q = e'^^' , k = -ie^'"'' , 6 > and a G R , 
the kernel of constant r-matrix is given by [1] 



(100) 



V(ai)(^(cr2)v?(cr3) 



'P{(^'l)vW2)'P{<^'?.) 
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where ^p{z) is the "non-compact quantum dilogarithm" [7] defined by 

(lUl) ip[z) - exp^^ J^_^.^smh{wb)smh{w/b) w ) ' ^(z + i6±V2) "1 + ^ 
Crossing parameter rj in (jlOOp is given by 

(102) , = ^. 

In Fock space {F~^) and anti-Fock space {F^) representations representations 

(103) q = e-^, fc = g"+V2 ^ n = 0, 1,2,3... (F+) or n = -1, -2, -3, . . . (F") , 
the matrix elements of constant r-matrix are given by a similar formula [1,2,21], 



3 



(104) 



where 



i=l 



^ 27ri J z"2+i (-g+"i-"3z;g2)^(_g-"i+"3^;g2) 



(105) r, = M2f)f if n = 0,l,2,3...(F^ 



or 



„n'{n'+l)(fj2. „2\ 

W ^.>-^ ;..(„.„(^;^)^: « n=-l,-2,-3,-4...(F-) 

The clockwise integration loop in (jl04|) circles all poles of the integrand but not includes 
z = 0. Pochhammer's symbols and Euler's quantum dilogarithm are defined by 

(107) {z; = (l - z){l - q'z) •••(!- g^C-D^) , (""/'^'f ^^ =l + z/q. 

{-qz;q^)oo 

Matrix p04|) has the block-diagonal structure in 

(108) Fl' F^^ <S) F^' , ei = ±, 

and thus it defines eight different i?-matrices. Classical limits of Fock and anti-Fock space 
representations (|103p provide 

(109) 0<k<l for F+ and Kk for F~ . 

Pre-factors and integrands in both pUOp and p04p have identical difference properties 
(leftmost relations in (jlOip and (|107p ). the main difference is that there is a non-compact set 
of poles in modular integrand and there is a compact set of poles in Fock space integrand. 
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Advantage of the special case ui = U2 = = 1 is that the constant r is the symmetric 
root of unity, 

(110) rj^s = 1 , r-|23 = SiS2S3T-l23('5l'52S3)"^ , 

where s = 1 for modular representation and Fock representation and 

(111) s = (-)" for F- . 

Factor s takes into account the anti-unitarity (a^)^ = — of anti-Fock representations. For 
instance, matrix 

(112) r'las = i-rr,23 in F+ ^ ^ F+ 
is Hermitian. 

Spectral parameters in (j96[) are given by "external field" factors. All cases below correspond 
to spectral parameters (132142 |49|55p with positive signs ei = e2 = £3 = 1. 

5.1. Regime 1. i?-matrix of ()96l) in modular representation < (I99p and spectral param- 
eters of Regime 1 is given by 

(113) i^l23 = e-2'''^2^Vi23e2'?<^i'^i+2'''^3'T3 _ 

Kernel of (jll3p is real and in the vicinity of equilibrium point (|35p it is positive with the 
asymptotic 

(114) (<7i(T2fT3|i?i23|o-i'720-3) ~ exp ^-^^ as 6 ^ and - \e'^^'^' ^ ki, 

where free energy Fq as function of (pi. .3 is given by (j39p . Presumably, partition function per 
site for cubic lattice and i?-matrix ()113p in physical regime < (3j < tt is 

(115) z = exp (^-^^0 
for arbitrary 77 > (I102p . 

5.2. Regime 2. i?-matrix of (j96p in modular representation fc^ < (|99p and spectral param- 
eters of Regime 2 is given by 

(116) i^l23 = ^-ie2''?<^2^Vi23e-2'''^i'^i-2'''^3^3 ^ 

where ^> is a unitary constant multiplier. In the vicinity — ie'^''""* — > fcj of ground state (j35p for 
Regime 2 the kernel of constant r-matrix oscillates. i?-matrix (jll6p is unitary and therefore it 
is the building block for a Heisenberg evolution operator. However, a spectral equation for the 
evolution operator is not yet known and we can't rigorously deduce a relation between spectra 
of quantum field theory and solitons and dispersion relation (I83p of classical field theory. 
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5.3. Regime 3. A self-consistent quantum field theory for Fock space representations corre- 
sponding to spectral parameters (I49|) and (jSOh is defined by 



(117) Rl23 = ^j-le'(^-'^2)n2^^2ggi0ini+i</,3n3 p- ^ p+ ^ p- _ 

Constant r-matrix in F.^ oscillates, operator ()117p is the unitary one for unitary 
constant multiplier g. A spectral equation for Heisenberg evolution operator is not known 
either except for a special 1 + 1 dimensional case and small occupation numbers [19]. In 
the same way as for Regime 2 we can not deduce rigorously relations between quantum 
spectra and classical dispersion relation. Note however an interesting feature of Regime 3: 
a self-consistency prescribes a correspondence between signatures of spectral parameters 
and a choice of representation F~^ with < /c < 1 or F~ with k > 1. Presumably, this 
correspondence provides a proper physical interpretation of spectra of evolution operators 
(see [19] for 1 + 1 dimensional case). A choice of constant g both in Regimes 2 and 3 is also 
a subject of proper physical interpretation. 

5.4. Regime 4. Curiously, the field-theoretical Regime 4 has no quantum field-theoretical 
counterpart; it corresponds to divergent statistical mechanics. 
i?-matrix for Regime 4 (I55p is given by 



(118) Ri23 = e'^2n2^'^^^g-0ini-02n2 p+ ^ p- ^ p+ 

where r' is defined by (jll2p and representation F^ F2 -^3^ is chosen in accordance with 
(I56p . Matrix elements of (jllSp are strictly positive and diverge as 

(119) ^ ^"ma-nma+nana+lower terms ni,n3^00, ^ -OO . 

A well-defined statistical mechanical lattice theory should involve a compensation of quadratic 
exponential asymptotic. It is possible via certain non-linear boundary conditions preserving 
the integr ability and involving extra three temperature- like parameters. 

In quasi-classical limit q = — > 1 diagonal matrix element of r' is given by 

(120) {n\r'\n) ~ Q<l>i^^-<f>2n2+<f>-inz {Yl\ as e ^ , 



^ e ^ 

where finite ki = g"' define by (j56p and (|57p the hyperbolic triangle with dihedral angles 9i 
and positive sides Vq is then given by (f59l) . For so defined 4>i the field factor in ()118|) 
compensates the pre-exponent in (I120p . However, this "classical equilibrium point" has no 
relation to a self-consistent quantum model. 

5.5. Regimes 2 and 3 as gauge field theories. Well defined quantum field theories in 
Regimes 2 and 3 involve Bose g-oscillators. However, the algebraic approach to the quan- 
tum tetrahedron equations allows one to introduce Fermi oscillators in addition to Bose 
ones [21]. All fermionic i?-matrices are even, they involve two fermions and one boson. Both 
Fermionic and bosonic ii-matrices are building blocks of Heisenberg evolution operator. One 
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can straightforwardly consider an evolution of simple test states with small total occupation 
numbers (in the Fock space representation for bosons). The evolution produces a set of Feyn- 
mann diagrams on constant time discrete surface (kagome lattice). In addition to simple 
propagation, fermionic i?-matrices are responsible for emissions of bosons, decay of boson 
into fermion pair and annihilation of fermion pair into a boson. Thus, the interpretation of 
quantum field theories as gauge field theories, where the bosons are gauge fields and fermions 
are matter field, is quite natural. Presumably, a proper choice of spectral parameters pro- 
vides also a gap between bosonic ground state and fermionic ground state, thus the spectral 
parameters are responsible in addition for fermionic mass. 

6. Discussion: Algebraic curves of higher genera 

Formulas (j68p for the soliton solution from the previous sections formally coincide with 
those for a general complex algebraic geometry (finite gap) solution: Pi,Qi are divisors on a 
genus g algebraic curve F^, £' is a prime form on it, / is related to a point on Jac(Fg), 0„ is 
a Theta-function: 

-■Pi /•P2 

:i2i) a 



n 



rfi r-f2 r±^3 

6(/(n)) , I{n) = z + ni / u - n2 I u + I u £ Jac(Fg), 

JQi JQ2 JQs 



where cj is a vector of Abel' holomorphic differentials and z is an arbitrary point on the 
Jacobian. Any of three-terms relation ([5]) is just the Fay identity [16]. Expression for I{n) 
in (jl2ip corresponds to a special case of homogeneous divisors. Divisors Pi,Qi are not free, 
they are divisors of three meromorphic functions, 

(122) N I u = mod(7r,7rO), i = 1,2,3, 



where is a size of cubic lattice, and is a period matrix; equations ()122p provide the 
periodical boundary conditions. 

Soliton solution is not just a straightforward trigonometric limit of algebraic geometry one 
since conditions ()122p are relaxed. Pi and Qi in general complex soliton solution are free. 

For the discrete time evolution system, the initial data of Cauchy problem define uniquely 
the algebraic curve [13] and thus selects the cases of soliton or finite gap dynamics. 

For general finite gap dynamics the spectral parameters Ui in ([5]) are not uniquely defined. 
The reason is that equations of motion ([5]) have gauge invariance. The gauge transformation 

(123) af^^ - 
such that 

(124) 6,Ti+e2 = 6,nC3,Ti 

is equivalent to the transformation of spectral parameters 

(125) Ul Ui- , U2 ^2— , Us U3- 
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Existence of homogeneous solution of equations of motion Q fixes the gauge group element 
and provides thus the definition of Uj. Otherwise, parameters Uj are irrelevant. 

Based on the principles of quantum-classical correspondence, one can conclude that the 
canonical quantization of g-oscillators (j94p and the choice of the Hilbert space as the product 
of local irreducible representations of g-oscillators corresponds to the choice of soliton sector 
on classical equations of motion. In particular, the condition of polynomial structure of 
Q-operators for a nested Bethe Ansatz for Fock space representations literally corresponds 
to factorization of spectral curve. The finite gap sector must correspond thus to another 
quantization scheme - a finite-gap quantization. 

Acknowledgements. I am grateful to V. Bazhanov, R. Kashaev, V. Mangazeev and P. 
Vassiliou for valuable discussions and fruitful collaboration. 
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